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A.  STATEMENT  OF  THE  PROBLEM  STUDIED 

Methods  for  representing  multi-dimensional  objects,  such  as  functions  of  several  vari¬ 
ables  and,  more  generally,  (hyper-)surfaces  is  the  main  objective.  One  goal  of  such  represen¬ 
tation,  whether  approximate  or  exact,  is  the  efficient  evaluation  of  the  object:  Multivariate 
Polynomial  Interpolation  as  well  as  Scattered  Data  Approximation  both  fall  into  .this  cat¬ 
egory.  Another  goal  is  a  representation  that  allows  one  to  identify  and  access  easily  and 
simultaneously  relevant  aspects  of  the  object.  The  topics  of  wavelets  and  Wey I- Heisenberg 
systems  belong  to  this  second  category. 

As  to  multivariate  polynomial  interpolation,  we  hope  ultimately  to  duplicate  the  suc¬ 
cess  of  univariate  polynomial  interpolation  as  a  basic  tool  in  scientific  computing.  Current 
focus  is  the  derivation  of  error  formulae,  using  a  multivariate  divided  difference  just  devel¬ 
oped. 

Scattered  data  approximation  is  to  be  accomplished  by  extending  well-known  and 
efficient  techniques  for  fitting  to  data  on  uniform  meshes. 

Concerning  Wavelets  and  Weyl-Heisenberg  systems,  the  goal  is  to  make  new  imoads 
in  these  important  areas  by  applying  tools  and  techniques  from  Approximation  Theory, 
particularly  those  developed  during  studies  of  Shift-Invariant  Spaces,  a  centerpiece  of  our 
previous  research. 

B.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

This  brief  outline  relies  on  the  fact  that  more  details  can  be  found  in  the  (semi-)annual 
reports  submitted  during  the  grant  period,  and,  if  need  be,  in  the  manuscripts  filed  with 
ARO  as  required.  All  bracketed  mimbers  refer  to  items  in  the  list  of  publications  given  in 
C  and  in  the  Bibliography. 

A  new  theory  for  wavelet  frames  was  developed  in  [13].  The  theory  is  based  on  tools 
and  observations  that  were  developed  in  the  context  of  shift-invariant  spaces^  and  in  partic¬ 
ular  the  ‘fiberization’  of  the  basic  operators  associated  with  a  shift-invariant  system,  [20], 
[22]  (as  part  of  work  on  an  earlier  grant).  The  new  theory  of  wavelet  frames  is  based  on 
the  interplay  between  the  given  wavelet  system,  and  a  new,  closely  related,  shift-invariant 
system:  a  quasi-wavelet  system.  The  theory  led  to  new  algorithms  for  constructing  wavelet 
frames  from  Multiresolution  Analysis.  The  highlight  of  these  algorithms  is  that  there  is 
essentially  no  restriction  in  them  on  the  scaling  function.  As  a  result,  compactly  supported 
spline  tight  frames  were  constructed  (in  one  and  in  several  variables)  [16].  The  new  frames 
may  be  useful  for  various  wavelet  applications,  most  notably  in  the  areas  of  feature  de¬ 
tection  and  noise  removal.  A  study  of  the  application  of  the  aforementioned  theory  was 
initiated  during  the  duration  of  this  grant. 

Several  other  studies  in  the  interface  between  Approximation  Theory  and  Wavelets 
were  successfully  completed.  Of  note  is  the  theory  of  approximation  orders  of  finitely 
generated  shift-invariant  (FSI)  spaces  (cf.  [5]).  The  theory  unravels  some  of  the  previous 
‘mysteries’  in  the  studies  of  wavelets  (such  as  the  so-called  ‘sum-rules’),  and  was  applied  in 
the  search,  [10],  of  new  factorization  techniques  of  wavelet  masks  in  one  dimension.  The 
connection  between  the  approximation  orders  of  wavelets  and  their  smoothness  was  also 
studied,  [11]. 


A  better  understanding  of  the  error  formula  in  the  proposers’  least  polynomial  in- 
terpolant  was  reached  for  special  cases  (see  [4]),  and  a  matlab  code  for  the  construction 
and  evaluation  of  the  least  interpolant  was  extended  to  arbitrary  dimensions  and  made 
available  on  the  web  (see  http://www.cs.wisc.edu/~deboor/multiiiit/m_files.html). 

Substantial  progress  has  been  made  in  the  area  of  approximation  to  scattered  data. 
First,  a  new  method  was  developed,  [21],  that  allows  the  conversion  of  the  (readily  avail¬ 
able,  and  very  well-understood)  approximation  schemes  on  uniform  grids  to  schemes  suit¬ 
able  for  scattered  data  approximation.  A  concrete  algorithm  based  on  the  general  conver¬ 
sion  method  was  developed  in  the  Ph.D.  thesis  of  Jungho  Yoon.  The  algorithrn  deals 
with  approximation  to  scattered  data  on  bounded  planar  domains.  Numerical  experiments 
with  a  Matlab  code  developed  for  that  purpose  showed  that  the  algorithm  performs  at  least 
as  well  as,  and  in  many  instances  better  than,  the  state-of-the-art  methods  for  scattered 
data  approximation. 

Some  progress  was  made  on  the  problem  of  proper  parametrization  of  curves  for  use 
in  the  construction  of  surfaces  by  blending,  in  that  Scott  Kersey  completed  a  thesis  on 
the  construction  of  a  smoothest  curve  that  fits,  within  given  tolerances,  to  a  given  point 
sequence  in  any  number  of  dimensions,  with  the  smoothness  measure  being  minimized  also 
with  respect  to  the  parametrization  used. 
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